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Abstract: 

We solve the problem of formulating Brownian motion in a relativistically covariant 
framework in 1 + 1 and 3 + 1 dimensions. We obtain covariant Fokker-Planck equations 
with (for the isotropic case) a differential operator of invariant d'Alembert form. Treating 
the spacelike and timelike fluctuations separately, we show that it is essential to take into 
account the analytic continuation of "unphysical" fluctuations in order to achieve these 
results. 
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Introduction 

Nelson 1 , in 1966, constructed the Schrodinger equation from an analysis of Brownian 
motion by identifying the forward and backward average velocities of a Brownian particle 
with the real and imaginary parts of a wave function. He pointed out that the basic process 
involved is defined nonrelativistically, and can be used if relativistic effects can be safely 
neglected. The development of a relativistically covariant formulation of Brownian motion 
could therefore provide some insight into the structure of a relativistic quantum theory. 

In recent years, Brownian motion has been applied as a mechanism for the description 
of irreversible evolution of quantum mechnical states, for example, collapse of wave func- 
tions as result of measurement. One finds that the wave function with a stochastic term 
added to the Schrodinger equation (represented in the Hilbert space of states 2 or in their 
projective representation 3 ) evolves to a mixed state of alternative outcomes of a measure- 
ment with the same probability distribution as given by the calculation of a priori (Born) 
probabilities in the framework of the standard quantum theory 2,3 . One of the motivations 
for the work of Adler and Horwitz 3 was based on the existence of a statistical mechanics, 
developed by Adler 4 and Adler and Millard 5 for the description of the equilibrium state 
of a general class of quantum fields; in this theory, the thermal average of commutator 
expressions take on the values implied by the familiar complex canonical quantum theory 
and quantum field theory. The Brownian motion corrections to the standard theories may 
be thought of as arising from the fluctuations around the equilibrium state. The develop- 
ment of a relativistically covariant theory of Brownian motion could make an extension of 
these ideas to relativistic quantum theory and quantum field theory more accessible. 

The program of stochastic quantization of Parisi and Wu 6 assumes the existence of 
relativistic wave equations, and applies a statistical approach closely related to path in- 
tegrals. The Wiener distribution of nonrelativistic Brownian motion, as is well known, is 
associated with (imaginary time) path integral formulations, and one would expect that 
covariant Brownian motion would similarly be associated with Parisi- Wu stochastic quan- 
tization. 

Nelson 1 has pointed out that the formulation of his stochastic mechanics in the context 
of general relativity is an important open question. The Riemannian metric spaces one 
can achieve, in principle, which arise due to nontrivial correlations between fluctuations in 
spacetime directions, could, in the framework of a covariant theory of Brownian motion, 
lead to spacetime pseudo- Riemannian metrics in the structure of diffusion and Schrodinger 
equations. 

In this paper we shall study the structure a covariant theory of Brownian motion. 

We first point out some of the obvious difficulties in reaching a covariant theory of 
Brownian motion, and indicate the directions we have chosen to solve these problems. 

Brownian motion, thought of as a series of "jumps" of a particle along its path, neces- 
sarily involves an ordered sequence. In the nonrelativistic theory, this ordering is naturally 
provided by the Newtonian time parameter. In a relativistic framework, the Einstein time 
t does not provide a suitable parameter. If we contemplate jumps in spacetime, to ac- 
comodate a covariant formulation, a possible spacelike interval between two jumps may 
appear in two orderings in different Lorentz frames. We therefore adopt the invariant 
parameter r introduced by Stueckelberg 7 in his construction of a relativistically covariant 
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classical and quantum dynamics. For the many body theory, Piron and Horwitz 8 postu- 
lated that this parameter is universal, as the Newtonian time in the nonrelativistic theory, 
and drives the classical particle trajectories Xi^(r) (worldlines labelled i = 1, 2, 3, ....AT) 
through equations of motion, and the evolution of the wave function in the Hilbert space 
L 2 (R 4N ) : ip T ({xi^}) through the Stueckelberg-Schrodinger equation (the differential form 
of the action of a one-parameter unitary group with parameter r). 

A second fundamental difficulty in formulating a covariant theory of Brownian mo- 
tion lies in the form of the correlation function of the random variables of spacetime. 
The straightforward generalization of the usual Brownian correlation property to special 
relativity, i.e., 



contains the serious problem that < dwo(r)dwo(T) >< 0, which is impossible. Brownian 
motion in spacetime, however, should be a generalization of the nonrelativistic problem, 
constructed by observing the nonrelativistic process from a moving frame according to the 
transformation laws of special relativity. Hence, as a first step, the process taking place 
in space in the nonrelativistic theory should be replaced by a spacetime process in which 
the Brownian jumps are spacelike. The pure time (negative) self-correlation does therefore 
not occur. In order to meet this requirement, we shall use a coordinatization in terms of 
generalized polar coordinates which assure that all jumps are spacelike. In this case, one 

would expect a distribution function of the form e «^ ; where [i is the invariant spacelike 
interval of the jump, and a is some constant. As we shall see, a Brownian motion based on 
purely spacelike jumps does not yield the correct form for an invariant diffusion process. 
We must therefore consider the possibility as well that, in the framework of relativistic 
dynamics, there are timelike jumps. In a frame in which the timelike jumps are pure 
time, the construction of the Gaussian distribution from the central limit theorem can 
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again be applied. The distribution would be expected to be of the form e - ^, where a 
is the timelike interval of these jumps, and b is some constant. By suitably weighting the 
occurrence of the spacelike process (which we take for our main discussion to be "physical" , 
since its nonrelativistic limit coincides with the usual Brownian motion) and an analytic 
continuation of the timelike process, we show that one indeed obtains a Lorentz invariant 
Fokker-Planck equation in which the d'Alembert operator appears in place of the Laplace 
operator of the 3D Fokker-Planck equation. One may, alternatively, consider the timelike 
process as "physical" and analytically continue the spacelike ( "unphysical" ) process to 
achieve a d'Alembert operator with opposite sign. 

2. Brownian motion in 1+1 dimensions 

We consider a Brownian path in 1 + 1 dimensions of the form 



We start by considering the second order term in the series expansion of a function of 
position of the particle on the world line, /(x M (r) + Ax^), involving the operator 





(2.1) 



O = Ax^Ax 



d d 



(2.2) 



dx^ dx v 
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We have remarked that one of the difficulties in describing Brownian motion in spacetime 
is the possible occurrence of a negative value for the second moment of some component 
of the Lorentz four vector random variable. If the Brownian jump is timelike, or spacelike, 
however, the components of the four vector are not independent, but must satisfy the 
timelike or spacelike constraint. Such constraints can be realized by using parametrizations 
for the jumps in which they are restricted geometrically to be timelike or spacelike. We 
now separate the jumps into spacelike jumps and timelike jumps accordingly, i.e., for the 
spacelike jumps, 

Ax — \x cosh a At = \i sinh a (2-3) 

and for the timelike jumps, 

Ax = a sinh a At = a cosh a (2-3') 

Here we assumed that the two sectors have the same distribution on the hyperbolic variable. 
We first look for the effects of a particle experiencing spacelike jumps only. In that case 
the operator O takes the following form: 



q2 q2 q2 

O = fj, 2 [cosh a— 7T + 2 sinh a cosh ol——— + sinh (2-4) 



If the particle going under timelike jumps only we find the operator O takes the following 
form: 

d 2 d 2 d 2 

O = a 2 [sinh 2 a-—r + 2 sinh a cosh a - - + cosh 2 a—A (2.5) 
ox 2 oxot at 1 

In order to obtain the relativistically invariant d'Alembert diffusion operator, the 
expression obtained in the timelike region must be subtracted from the expression for 
the spacelike region, and furthermore, the amplitudes must be identified. In the physical 
timelike region, the coefficient a 2 is, of course, positive, and using the law of large numbers 
on the random distribution, one obtains a Gaussian distribution analogous to that of the 
spacelike case. 

We see, however, that we can obtain the d'Alembert operator only by considering the 
analytic continuation of the timelike process to the spacelike domain. This procedure is 
analogous to the effect, well-known in relativistic quantum scattering theory, of a phys- 
ical process in the crossed (t) channel on the observed process in the direct (s) channel. 
Although we are dealing with an apparently classical process, as Nelson has shown, the 
Brownian motion problem gives rise to a Schrodinger equation, and therefore contains 
properties of the quantum theory. We thus see the remarkable fact that one must take 
into account the physical effect of the analytic continuation of processes occurring in a 
non-physical, in this case timelike, domain, on the total observed behavior of the system. 

In the non-stochastic case, Einstein's relativity identifies Ax/ At with p/E, where p 
and E are the components of the energy-momentum four- vector of the particle. If we make 
an analogous identification, assigning these variables as properties of the fluctuations, then 

a 2 = (At) 2 - (Ax) 2 oc (E 2 - p 2 ), (2.7) 
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defining a stochastic mass squared associated with the Brownian particle. If the relation 
between E and p becomes spacelike, the notion of stochastic mass can be retained undere 
the transformation to an imaginary representation E — > iE' and p — > ip' , for E',p' real§, 
so the relation p' / E' > 1 remains, but p' 2 — E' 2 > 0. We assume sufficient symmetry in 
the spacelike and timelike distributions so that the absolute values of < a 2 > and < /j 2 > 
are equal. The preservation of the mean magnitude of the interval reflects the conservation 
of a mass-like property which remains, as an intrinsic property of the particle, for both 
spacelike and timelike jumps. 

With these assumptions, the cross-term in hyperbolic functions cancels in the sum, 
which now takes the form 

<°> =< " 2 > l& - & (2 - 8 » 

Taking into account the drift term in (2.1), one then finds the relativistic Fokker-Planck 
equation 

where d/dx^ operates on both and D. 

We see that the procedure we have followed permits us to construct the Lorentz invari- 
ant d'Alembertian operator, as required for obtaining a relativistically covariant diffusion 
equation. Furthermore, since the expectation of sinh 2 a, cosh 2 a could be infinite (e.g., 
for a uniform distribution on a), the result we obtain in this way constitutes an effective 
regularization. 

3. Brownian motion in 3 + 1 dimensions 

In the 3 + 1 case, we again separate the jumps into timelike and spacelike types. The 
spacelike jumps may be parametrized, in a given frame, by 

At =\i sinh a 

Ax =u cosh a cos 6 sin ft 

(3-1) 

Ay =n cosh a sin <fi sin i? 
Az =ji cosh a cos ■& 

We assume the four variables n, a, are independent random variables. In addition 
we demand in this frame that $ and <p are uniformly distributed in their ranges (0, 7r) and 
(0, 2ir), respectively. In this case, we may average over the trigonometric angles, i.e., $ and 
(f> and find that: 

H 2 

< Ax 2 > <t>i # =< Ay 2 >^=< Az 2 >4>,&= Y cosh2a ^3 2 ) 

< At 2 >^ i ^= / u 2 sinh 2 a 



§ This transformation is similar to the analytic continuation p — > ip' in nonrelativistic 
tunneling, for which the particle appears as an instanton. 
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We may obtain the averages over the trigonometric angles of the timelike jumps by replac- 
ing everywhere in Eq.(3.2) 



2 2 

cosh a <-> sinh a 

2 2 
[l — > (T 



to obtain 



< Ax >^ =< Ay ><^=< Az ><^= y sinh a 

< At 2 >^ ; ^=cr 2 cosh 2 a, 

where a is a real random variable, the invariant timelike interval. Assuming, as in the 1 + 1 
case, that the likelihood of the jumps being in either the spacelike or (virtual) timelike 
phases are equal, and making an analytic continuation for which a 2 — > —A 2 , the total 
average of the operator O, including the contributions of the remaining degrees of freedom 
/z, A and a is 

d 2 

<0 >=(< a 2 >< sinh 2 a > - < A 2 >< cosh 2 a >)tt^+ 

dt (3.4) 

- (< fi 2 >< cosh 2 a > - < A 2 >< sinh 2 a >) A 

If we now insist that the operator < O > be invariant under Lorentz transformations (i.e. 
the d'Alembertian) we impose the condition 

< fi 2 >< sinh 2 a > - < A 2 > < cosh 2 a >= 

l/o o o 9 \ (3-5) 

- - (< /i >< cosh a > - < A >< sinh a >) 

Using the fact that < cosh 2 a > — < sinh 2 a >= 1, and defining 7 =< sinh 2 a >, we find 
that 

< * 2 >= < " 2 > (3 ' 6) 

The Fokker-Planck equation then takes on the same form as in the 1 + 1 case, i.e., the 
form (2.9). We remark that for the 1 + 1 case, one finds in the corresponding expression 
that the 3 in the denominator is replaced by unity, and the coefficients 4 are replaced by 
2; in this case the requirement reduces to < /j 2 >=< A 2 > and there is no 7 dependence. 

We see that in the limit of a uniform distribution in a, for which 7 — > 00, 
< A 2 >^< fj 2 >. In this case, the relativistic generalization of a nonrelativistic Gaussian 

distribution of the form e <« becomes of the form e 37, which is Lorentz invariant. As in 
the 1 + 1 case, the result (3.6) corresponds to a regularization. 

4. Conclusions and Discussion 

We have constructed a relativistic generalization of Brownian motion, using an invari- 
ant world-time to order the Brownian fluctuations, and separated consideration of spacelike 
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and timelike jumps to avoid the problems of negative second moments which might other- 
wise follow from the Minkowski signature. Associating the Brownian fluctuations with an 
underlying dynamical process, one may think of 7 in the 3 + 1 case as an order parameter, 
where the distribution function (over a), associated with the velocities, is determined by 
the temperature of the underlying dynamical system (as we have remarked, the result for 
the 1 + 1 case is independent of the distribution on the hyperbolic variable). 

At equilbrium, where dD/dr = 0, the resulting diffusion equation turns into a classical 
wave equation which, in the absence of a drift term K^, is the wave equation for a massless 
field. An exponentially decreasing distribution in r of the form exp — kt would correspond 
to a Klein-Gordon equation for a particle in a tachyonic state (mass squared —k). We 
have considered the spacelike jumps as "physical" since they result in the usual Brownian 
motion in the nonrelativistic limit. If the timelike jumps were considered as "physical" , 
one would analytically continue the "unphysical" spacelike process. The resulting diffusion 
equation would have the opposite sign for the d'Alembert operator, and an exponentially 
decreasing distribution would then result in a Klein-Gordon equation in a timelike particle 
state. 

Nelson 1 has shown that non-relativistic Brownian motion can be associated with a 
Schrodinger equation. Equipped with the procedures we presented here, constructing rel- 
ativistic Brownian motion, Nelson's methods can be generalized. One then can construct 
relativistic equations of Schrodinger (Schrddinger-Stueckelberg) type. The eigenvalue equa- 
tions for these relativistic forms are also Klein-Gordon type equations. Moreover one can 
also generalize the case where the fluctuations are not correlated in different directions 
into the case where correlations exist, as discussed by Nelson 1 for three dimensional Rie- 
mannian spaces. In this case the resulting equation is a quantum equation in a curved 
Riemannian spacetime; as pointed out in ref.10, the eikonal approximation to the solu- 
tions of such an equation contains the geodesic motion of classical general relativity. The 
medium supporting the Brownian motion may be identified with an "ether" 11 (Nelson 1 
has remarked that the self-interaction of charged particles might provide a mechanism for 
the Brownian motion) for which the problem of local Lorentz symmetry is solved. This 
generalization of Nelson's method will be discussed elsewhere. 

Acknowledgements 

We would like to thank Linda Reichl, W.C. Schieve and Sungyan Kim at the Ilya 
Prigogine Center for Statistical Mechanics and Complex Systems at the University of 
Texas at Austin for helpful and stimulating discussions. 

References 

1. Edward Nelson, Dynamical Theories of Brownian Motion, Princeton University Press, 
Princeton (1967); Edward Nelson, Quantum Fluctuations, Princeton University Press 
Princeton (1985). See also, Ph. Blanchard, Ph. Combe and W. Zheng, Mathematical 
and Physical Aspects of Stochastic Mechanics, Springer- Verlag, Heidelberg (1987), for 
further helpful discussion. 

2. For example, G.C. Ghirardi, P. Pearle and A. Rimini, Phys. Rev. A42, 78 (1990). 

3. L.P. Hughston, Proc. Roy. Soc. London A452, 953 (1996), and references cited there; 
S.L. Adler and L.P. Horwitz, Jour. Math. Phys. 41, 2485 (2000). 



7 



4. S.L. Adler, Nuc. Phys. B415, 195 (1994); S.L. Adler, Quaternionic Quantum Median 
ics and Quantum Fields, Oxford Univ. Press, N.Y. (1995); S.L. Adler, [hep-th/0206T20 

5. S.L. Adler and A.C. Millard, Nuc. Phys. B473, 199 (1996). 

6. G. Parisi, Y. Wu, Sci.Sin. 24 483 ; Mikio Namiki Stochastic Quantization , (Springer- 
Verlag, Heidelberg, 1992). 

7. E.C.G. Stueckelberg, Helv. Phys. Acta 14, 322 (1941); 14, 588 (1941). 

8. L.P Horwitz and C. Piron, Helv. Phys. Acta 46, 316 (1973). 

9. See, for example, S.S. Schweber, An Introduction to Relativistic Quantum Field The- 
ory, Harper and Row, N.Y. (1961). 

10. L.P.Horwitz , O.Oron, |hep-ph/02050i~8 

11. S. Liberati, S. Sonego and M. Visser, Ann. Phys. 167 298 (2002). 



8 



